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A  GENERALIZATION  OF  THE  K-V  EQUILIBRIUM  FOR  AN  INTENSE 
RELATIVISTIC  BEAM  PROPAGATING  IN  A  SOLENOID 


1.  nmtoDucnm 

The  Vlasov  equilibrium  for  a  charged  particle  beam  in  a  uniform 
axial  magnetic  field  with  distribution  function 

f(r,p)  .  nQ/2xrQ«  i(H-«L-YQnc^)  8(P-fl^Yo«c) 


is  examined,  vhere  n,.,  Yn>  ^  constants.  An  expansion  for 

^  2  2  2  2  2 
which  the  dimensionless  quantities  4)tb  e  Oq/y^bc  ,  «tb/c,  and  1-^  >  I/Yq 

are  small  is  taken  idiere  b  is  an  order-one  scale  length.  In  the  limit 

that  Tq  '*  **  nnd  «  0,  the  K-V  distribution  is  recovered  {1,2,  pp.  588- 

594].  The  four  constants  in  the  distribution  function,  n^,  Yq>  n  and 

8^,  are  related  to  the  more  natural  specification  of  the  beam  by  its 

total  current,  total  energy,  canonical  angular  momentum,  and 

perpendicular  oiittance,  e^,  respectively.  As  an  example  of  a  case 

vhere  the  corrections  to  the  K-V  distribution  are  important,  ve  discuss 

the  beam  proposed  for  the  Spiral  Line  Induction  Accelerator  (SLIA) 

Proof-of-Concept  Experiment  (POCE)  [3]. 
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2.  CHOICE  OF  VLASOV  EQUIUBRIUH 


Consider  an  equilibrium  which  is  independent  of  z,  the  propagation 
direction,  and  9.  The  only  fields  are  E  (r)  «  -d*/dr,  B_(r)  «  1/r 
drA^dr,  and  Bg(r)  >  -dA^/dr,  where  #  is  the  scalar  potential  and  A^  and 
A^  are  components  of  the  vector  potential.  The  constants  of  motion  are 

2 

H  -  ymc  +  Zei,  the  energy, 

L  «  *  ZeAg/c),  the  canonical  angular  momentum,  and 

P  >  p  -f  ZeA  /c,  the  canonical  axial  momentum, 

[2,  pp.  84-88].  Ve  will  retain  the  charge  number  Z  which  is  -1  for 
electrons.  The  situation  we  are  most  interested  in  is  the  description 
of  a  bean  traveling  in  an  external,  uniform  solenoidal  field,  B^  -  B^ 
for  which  AQ***(r)  » 

Any  distribution  of  the  form  f(r,p)  «  fQ(H,L,P),  where  fQ  is  a 
given  function,  is  an  equilibrium  solution  of  the  Vlasov  equation  [2]. 

As  will  be  shown,  the  fora  from  idtich  the  R-V  can  be  derived  is  given  by 

f(r,p)  -  nQ/2nYQn  8(H-«L-YQnc^)  S(P-PjjYQmc)  (1) 

where  nQ,«»,YQ  and  ^  are  constants. 

2  2 

Denote  the  dimensionless  potentials  4  >  e#/nc  ,  a  -  eA  /me  and 

2  z  z 

»e  -  • 

Consider  an  integrable  function  g(r,p^,Pg,P2).  Define  <g(r,p)>  a 
/g(r,p)f(r,p)d^p,  then  it  can  be  shown  that 

<g>  -  dx(g(r  ,ncX  cos(T),acXsin(T)  +  ncv,nc(^YQ-Za2))  + 

in/2 

g(r,-ncX'cos(T),ncXsin(x)  ♦  mev,  mc(P^YQ-Za^)))<U  +  icsin(x))  (2) 
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vhere 


The  assuBption  that  |vr|  «  jur/cl  <  1  has  been  made. 

The  beam  radial  extent  is  determined  by  the  region  for  vhich  the 
argument  of  the  square  root  in  the  expression  for  X  is  non-negative. 

Ve  are  interested  in  parameter  regimes  vhich  allov  rapidly 
convergent  pover  series  in  the  perpendicular  momenta.  To  the  lowest 
orders  we  have  -  1  -  -  1/Yq^  +  order-one  constant  tines  (♦/yi 

a-/Y»  a  /y  or  u^r^)  and  ~  v^r^.  For  relativistic,  low  emittance 

beams  the  first  term  is  small  because  ^  ■  i  ond  Tq  »  1*  However,  an 
expansion  in  powers  of  the  perpendicular  momentum  is  also  valid  when  Yq 
-»  1  and  0  at  least  for  low  charge  density  beams  vhere  ♦/yq,  s^/yq 
and  a^/YQ  are  small. 
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Note  that  <g>  «  0  for  all  odd  functions  of  p^. 


3.  EXAMPLES  OF  MOMENTS  OP  THE  DISTRIBUTION 

1)  Particle  density;  n  »  <1>  gives 

n  "O 

«  “  (l-v2r2j3/2 

This  result  is  exact. 

For  o^O  and  f(r,p)  «  SCH-YqBC*^)  the  result  is  n^  =  Oq  (1  -  Z^/yq) 
[4,  2,  pp.  122-130]  and  the  scalar  potential  is  a  modified  Bessel 
function  out  to  the  bean  radius.  The  "nonrelativistic”  or  "rigid-rotor" 
case  is  recovered  if  the  u  terns  are  dropped  and  low  current  is  assuned 
so  that  Z4/yq  is  ignored.  This  results  in  the  flat-top  radial  density 
of  the  R-V  distribution  [2,  pp.  588-594].  The  beams  under  study  will 
not  be  flat-top,  but  they  vill  have  a  sharp  radial  cutoff. 

2)  Moments  of  the  nonenta;  All  integral  povers  of  the  momenta  can 
be  integrated  exactly.  The  ones  of  imediate  interest  are 

<Pg>/mc  ■  nQ  (Xk/2  +  pv) 

<P^/b^c^  «  Sq  (X^u/2  +  XvK  +  v^y) 

<p^/mV  -  Sq  (3X^ic/8  +  3X^vm/2  +  3Xv^ic/2  +  v^y)  (5) 

<pj>/n^c^  .  nQ  x'^y/2  -  nQ(l-v^r^)X^y/2 

<PQpJ>/n^c^  .  nQ  (Xx'^k/8  +  x'^vy/2)  -  nQ(l-u^r^)(X^K/8  +  X^yv/2). 

Moments  involving  y  result  in  elliptic  integrals. 


vrZa-  Z4  'j 

1  +  - 5 - .  (4) 

Yq  Yq  J 


t 


I 


1 


% 


’'A 
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Substituting  from  Eq.  (3),  ve  get 


^  .  _  ..2  2%.X/2  e  X  X  ..2«. 

<p^>/Yo«c  -  n^vr  (1-v  r  )  ;;  j-r  ^  w  ) 


2yo  (l-v‘r‘) 


urZag  Zf 


Consider  the  canonical  angular  momentum  density, 


r<pQ>  rZagO^  2 


Yomc  Yo«c 


r  HqU  +  rnp  — 


•t-  higher  order  terms. 


2  2 

In  a  solenoid  ra^  »  r  eB^/Zmc  if  the  self-fields  are  neglected.  In  the  case 

of  a  beam  created  in  a  diode  free  of  magnetic  field,  the  canonical  angular  mo- 

2 

mentum  is  sero.  This  is  assured,  to  lovest  order,  by  ta/c  «  v  «  -ZeB^/Zmc  Yq  * 
B^jCkG)/  3,4yq  cn‘^. 

2  2  2 
Denote  Pj^  •  ■*"  Pe* 


2  2  2  2 
<p^>/yjm-^c^  - 


Uq  r  virZaQ  Z^l  2  +  3v^r^ 

L  Yq  yJ  Z(l-vV) 

2  .  V^r^  f  1  L  ^ 

“  ■  g  '  2  I  rj  ‘  ■  r„ 

,  2.,  3  3  3  .  „  f  20  .  f, 

^  °  8(1  -  Yn  yj 


,  and 


^IrK-^lT 
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4.  EXPANSION  IN  THE  "BEAM"  APPROXIMATION 


In  order  to  conpute  an  equilibriun  configuration  for  this 
distribution,  ve  need  to  compute  several  aoeents.  Hence,  ve  nov  examine 
the  expansion  of  r  in  an  approximation  relevant  to  beams  of  interest. 
Define  Yj,  -  and  .  It  follows  that 

In  order  to  call  this  a  "beam”  ve  make  the  following  assumptions: 

^  ‘  ^  1  ' 

<Pr^  +  Pq^)  «  (Yjj«c)^  ,  and 


I  Ze(A,,A.,*)  /  Y^mc^  |  «  1  . 


The  last  assumption  can  result  from  a  "small”  current  or  a  "large"  y> 
The  relativistic  factor,  y,  is 

f  ,  P?  -  p|  ♦  (  '  - 


f  ,  p;  ‘  p;  ♦  ( 

■  * - 7? 

me 


2  2  2  2  2 

which,  to  first  order  in  (p^+  p^)  /  Yjj«  c  and  a^/Yjj  ,  can  be  written 


2  2 


Y  »  Yu 


/  ^  A 

1  Pr  ^  Pe 
^  i  2  2  ” 


where  ve  have  used  P  «  H^YqBC.  To  the  sane  order, 

^  “  ^b  I  lyVc^  *  Yk  ]  ’ 


-  ^b  k  2Y^m‘c‘-  Yb 

’i  p,  ..  iiliil  ^ 

r"  “  Y«c  T)  ■  2  2,2^2  ■  3  ’ 

V  ^b 


'e  Pe 


C  Y«C  Yb*C 


V  ^b 
2  „2 
,  Pr  ♦  Pe 

1 - *  ft  A"  ♦ 


2Yb»‘c‘ 
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ymc 


Yk«C 


1  - 


P  P  e 


f 

Yk 


and 


-2  2. 


— j-  -  ifb  ‘  -  — - rr 

■c  zc 


1  pI  +  (Pe-  rwYb*)^ 

^  ? - X5 - 

Yb  ■  c 


Note  that  by  keeping  the  self-field  corrections  in  y,  the  aean  azimuthal 
velocity  -  <v^/<l>  is  not  exactly  r«»  (referred  to  as  "rigid-rotor" 
equilibria). 


5.  CALCULATION  OF  SOURCE  TERMS  FOR  MAXWELL'S  EQUATIONS 

We  calculate  the  charge  and  current  densities  to  the  order 

2  2 

indicated  above:  ve  retain  tens  t  4/y>  uba^/Yt  u  b  ,  vbi/y, 

uba^/Y»  u^b^Sg/Y*  v^b^,  <1  -  Yb^/Y0^),  and  v)b(l  -  Yb^^YQ^)  where  b  is  a 
characteristic  length  of  the  saae  order  as  the  radial  diaension  for  the 
bean.  The  value  of  b  is  uniaportant,  as  it  does  not  appear  in  the  final 
equations. 


1.  Charge  density,  p  >  Ze<l>  : 

p  -  Zenp  (  1  ♦ 


2  2 

3v  r  vrZa, 


^0  ^0  ^ 


2.  Axial  current  density,  «  Ze<v^>  : 

.  ..2_2 


J,-  Zecno 


-fi. 

^12  2?)  2  I  y5  J  Yk 


f  '»z«e 

z* 

1  1  0 

“2 

''  ^0 

^0  ^ 

^  ^b  ^ 

>  Zccn« 


(10) 


(11) 


3.  Azimuthal  current  density,  >  Ze<v^: 


If  only  the  first  term  is  ki^pt,  ve  obtain  Jg  >  Zecn^vr,  ra^  becomes  a 

fourth-order  polynomial  in  r,  and  4  and  a  are  the  suns  of  fourth-order 

polynomials  and  modified  Bessel  functions  (inside  the  bean).  For  u  «  0, 

the  expression  for  J_  reduces  to  the  correct  expression  to  the  proper 

z 

order  in  the  expansion  (4,2,  pp.  122-130].  For  the  K-V  distribution, 
only  the  first  terms  in  the  expressions  for  p,  and  are  retained 
(2,  pp.  588-594]. 


Let  us  introduce  some  dimensionless  quantities.  The  dimensionless 

radial  coordinate  is  r/b.  The  dimensionless  particle  density  is  n,.  ■ 

2  2  2  2  ^ 
4nQnb  r^  •  nQb(cm)  /2.83ell  where  i®  the  classical  electron 

radius  (  a  uniform  1  kA  beam  of  1  cm  radius  has  a  particle  density  n  - 

5  •• 

6.63cl0  ca~  ).  The  dimensionless  angular  parameter  is  t»  ■  «9b/c.  The 

dimensionless  external  solenoid  field  is  Sq  ■  BQeb/nc  -  BQb/1.7kG.  The 

dimensionless  current  density  is  J  ■  J4xeb  /mc^  and  the  dimensionless 

beam  current  is  f  a  le/mc  ■  I(kA}/17.  Nov  the  equations  will  use 

dimensionless  quantities  and  drop  the  and  ve  will  use  »  instead  of 
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6.  POWER  SERIES  SOLUTION  FOR  THE  POTENTIALS 

We  now  use  the  power  series  expansions  of  the  potentials: 

♦  •  '  *2  “  ’^*0  "  Setting  the  potentials 

equal  to  0  at  the  origin  allot::  the  suns  to  run  from  n-1  to  n*<>>. 

We  also  define  coefficients  a„,  b  and  c  for  n^l  to  5  froe 

n  n  n 

Maxwell's  equations  for  the  potentials  inside  the  beam: 


1  d  ,  d*  ,  2 

?  3F  d?  -  -  *1  Ve  *  “s’"  ’ 

t  df  "  sf*  "  ’^z  •  •‘i  •»2^  ♦  Vz  *  Ve 

r  ''*0  -  -^0  -  ^  ^2*  *  Vz  *  Ve  ^  ^5^^^  ' 


where  p  is  the  dimensionless  charge  density  with  eng  replaced  by 
dimensionless  Sq,  v  replaced  by  dimensionless  »  *  ub  and  r  replaced  by 
r/b  in  Eq.  (10).  Similar  substitutions  are  made  in  the  previous 
expressions  for  J,  -»  and  J-  ->  Jq. 

Z  2  9  V 


The  following  equations  for  the  coefficients  result: 

4Sj  -  Sj  .  -Zuq  , 

•  "^0^0  ’ 

2  2  2 
aj  -  Z  u^Sj^/tq  -  Z  UqUU^/Yq  -  IZn^u  /2 

♦  bj  -  Z^n^tj/Yo  -  ZnQ^v2/2  , 

8u2  •  Cj  •  ~^0**  * 


4tj  -  bj 


f 
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240^  ■  *^2®!  ^3^1  ^5 

36S3  *  ^  ®3^2  ^  ®4^2 

and 


36t- 


^2®2  ^  ^3^2  *  ^4^*2 


-  2^nQ«Sj/rQ  -  z\q(0^u^/Yq  -  3ZnQ(i)^/2  , 

2  2 
>  Z  nQS2/YQ  -  Z  nQau2^YQ  > 

*  Z^nQt2/YQ  •  .  (14) 


Note  that  all  the  coefficients  are  deterained  explicitly  except 
which  will  be  found  later.  Coefficients  of  order  n  depend  on  those  of 
lover  order.  Specifically,  s^^  -  n^/y  ~  ^n-1’ 

and  u^  ~  n^w/Y  2'*‘***n  2^*  Therefore  the  power  series  converge 
rapidly  and  can  be  terainated  provided  Rq/yq  «  1.  Note  that  this 
condition,  which  is  essentially  the  requirement  for  ignoring  nonlinear 
terns  in  the  potentials,  is  different  than  the  condition  for  rapid 
convergence  in  expansions  in  perpendicular  momenta. 


From  now  on  we  consider  the  case  of  an  external  solenoid  of 
(dimensionless)  field  Bq  and  replace  u^  in  the  above  equations  with  u^ 
Bq/2  where  Uj^  from  now  on  refers  only  to  the  self-field  contribution. 

Ve  obtain 

2  Bgr  ^0**  4 
rag  -  uir  *  ^  r 


-I 


Z^nQ»  Z^nQ«#^Uj  Z^n^w^Bp  Zn^w^  g 

-56;- *  " -48;— "  J"  ' 


-ZnQ  2  "qWBq  4 

♦  -  —  r  ♦  - gj— -  r  ,  and 


(15) 


s  - 


2 

— '  r  ♦ 


64y/i 


32 
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We  will  need  the  following  expressions  to  the  proper  order 


«u:a 

♦ 


e 


-Znor‘ 


4Yn 


2^0  ’ 

4 

Z  n^w^r 


32yX 


and 


(16) 


z 

^0 


To  lowest  order  we  expect  Bq/Yq  "  w»  so  the  terms  in  4/y 
containing  r^  can  usually  be  ignored. 


7.  EXPRESSION  FOR  THE  BEAM  RADIUS 


The  beam  edge  radius  R  is  obtained  by  finding  the  domain  where  the 
argument  of  the  square  root  in  the  expression  for  X  in  Eq.  (3)(  is 
positive.  To  this  order  the  beam  is  confined  to  spatial  regions 
satisfying 


1  - - - j  +  - +  «  r  2  0 

Yq  ^0  ^0  ^0 


and 


we  assume  1  -  y^/Yq  >  0  so  that  the  beam  is  not  hollow.  The  outer 
cutoff  radius,  R,  of  the  beam  is  then  the  solution  of 


(17) 
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8.  CONSIDERATION  OF  CANONICAL  ANGULAR  MOMENTUM 


The  canonical  angular  BoiaentUB  is  L*(r)  -  <L>  >  <^(P0  *  Z^q/c)>. 
Expamding  this,  we  obtain 

L  /hQYQnc  ■  «4r^(3/2  -  Yjj^/2yq^  +  3«^ +  3ZcaraQ/YQ 

-  3Z4/Yo  +  ^Za^/Tg)  +  Zra^/Yo 

L  /npYginc  -  r^(  «(3/2  -  Yj,^/2yq^)  +  ZUj/yq  +  ZBq/2yq  )  + 

r^(  3«^  +  Z^Migt  3  -  Pq^J/Ayo 

+  3Z«^l  Uj  +  Bq/21/yq  ).  (18) 

For  a  non-iamersed  cathode,  as  is  the  case  for  the  SLIA  experiment, 

/2iirL*dr«0.  (19) 

Eq.  (19)  relates  the  parameter  m  to  the  external  field,  Yq>  8^*  and  u^. 
Other  diode  configurations  lead  to  different  values  for  the  beam's  total 
canonical  angular  momentum  [5,  2,  pp.  559-569]. 


9.  DETERMINATION  OF  THE  AXIAL  SELF-INDUCED  MAGNETIC  FIELD 

The  integration  of  Faraday's  lav  gives  the  self-induced 
(dimensionless)  axial  magnetic  field  as 

-  ligdr  -  ZnQ»«^/2  +  3ZnQ«V/8  +  Z^n^eR^/AYo  (Zop/A)  (20) 
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If  Yq.  and  Hq  are  chosen,  then  this  expression  for  Bq.  (17)  for 
the  bean  radius  R,  and  Eqs.  (18)  and  (19)  for  the  zero  canonical  angular 
■OBcntUB  give  three  (nonlinear)  equations  for  the  unknowns  Uj^,  R,  and  m, 
and  the  equilibriun  is  conpletely  specified. 

In  the  sinulations  to  be  presented  later,  the  bean  is  inside  a 

conducting  cylinder  of  radius  R„.  In  this  case,  the  relationship  for  u. 
4  2  2  2  * 

is  Uj^  «  -Uj  R  /R^  +  (1-R  /R^  )XjQdr/2.  This  expression  accounts  for 
the  axial  Magnetic  field  produced  by  the  iaage  currents. 

10.  SPECIFICATION  OF  THE  EQUILIBRIUM  IN  LABORATORY  PARAMETERS 

The  Vlasov  equilibriun  was  specified  in  terns  of  four  parameters, 
Hqi  Yq>  8^  snd  a  with  certain  restrictions  on  their  values  so  that  the 
expansion  is  valid.  Noraally,  however,  the  four  parameters  commonly 
chosen  are  the  bean  energy,  beam  current,  emittance  and  canonical 
angular  momentum. 

Ve  have  seen  that  conditions  on  the  canonical  angular  momentum  are 
closely  related  to  the  value  of  t».  Clearly,  bean  current  is  correlated 
with  Oq.  In  an  experiment  the  beam  electrons  are  created  with  an  energy 
determined  by  the  potential  on  the  anode  foil.  The  bean  usually 
propagates  inside  a  conducting  cylinder,  so  the  field  energy  must  be 
taken  from  the  beam.  Note  that  the  tube  boundary,  R^,  appears  here  if 
ve  wish  to  relate  the  diode  voltage  to. the  appropriate  value  of  yq* 

Ve  define  the  normalised  perpendicular  RMS  emittance,  e^,  in  the 
relativistic  setting,  throu|^ 

.  «r^»«Pp^  +  -  «rpj.»^  -  «rpg»^  (21) 

where 

«g(r,p^,PQ,p^»  a  Jgfd^p2lirdr  /  J£d^p2xrdr.  (22) 
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The  emittance  is  related  to  the  amount  of  the  beam's  energy  found  in 
"undirected”  transverse  motion  [6].  The  Vlasov  parameter  ^  is  a  measure 
of  hov  much  of  the  beam's  momentum  is  directed  axially  and  if  the  other 
parameters  are  fixed,  increasing  will  decrease  the  transverse  emittance. 

Figures  la-c  show  the  beam  current,  emittance  and  cutoff  radius,  R, 
as  functions  of  the  parameters  in  the  approximate  Vlasov  equilibrium.  Ve 
have  fixed  Yq  «  9.  Bq  takes  the  values  4.5  kG  (dashed),  5.5  kG  (solid), 
and  6.5  kG  (dotted).  is  set  to  .97,  .988  or  .993.  The  assumption  of 
zero  total  canonical  angular  momentum  is  made.  Note  that  the 
approximations  made  in  this  model  become  less  accurate  as  the  parameter  Hq 
increases  in  the  figures. 

11.  PROOF-OF-CONCEPT  EXPERIMENT  (POCE)  IN  THE  SPIRAL 
LINE  INDUCTION  ACCELERATOR  (SLIA) 

The  SLIA  is  a  proposed  compact  high-current  electron  beam  accelerator 
[3]  vhich  is  being  developed  and  studied  experimentally  at  Pulse  Sciences, 
Inc.  SLIA  uses  induction  cavities  to  accelerate  the  beam  and  a  solenoidal 
magnetic  field  for  transport  in  the  strais^t  sections.  Transport  around 
the  bends  is  accomplished  using  strong  focusing  stellarator  windings  (i.e., 
twisted  quadrupoles)  augmented  with  a  vertical  magnetic  field.  The  bean 
line  is  shown  schematically  in  Fig.  2.  The  transport  lines  are  isolated 
from  each  other  except  in  the  acceleration  region.  In  the  POCE  a  10  kA,  35 
ns  beam  will  be  accelerated  to  9.5  NeV  from  3.5  NeV  in  steps  of  1.5  MeV. 

The  motivation  for  this  paper  came  from  attempts  to  simulate  the 
POCE  experiment  [7].  Figure  3a  shows  a  simulation  of  a  10.7  kA,  r  -  9.3 
beam  with  a  normalized  RMS  emittance  of  .34  ca-rad  in  a  5.5  kG  solenoidal 
field.  The  beam  particles  are  loaded  with  a  K-V  distribution  in  the 
rotating  frame  with  the  correct  moments  for  a  matched  equilibrium. 
Nevertheless,  the  beam  radius  exhibits  growing  oscillations  at  the 
envelope  frequwicy.  In  this  example  the  (transverse)  emittance  actually 
decreases  slightly.  After  many  meters  the  beam  radius  oscillations  damp. 
This  behavior  indicates  that  the  K-V  equilibrium  is  not  a  good  kinetic 
equilibrium  for  these  beam  parameters. 
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The  next  figure  shows  the  results  of  a  simulation  with  the  beam 

initialized  with  the  equilibrium  developed  in  this  paper.  Specifically, 
12  3 

n^  «  1.4x10  per  cm  ,  Yq  >  9,  0^  «  0.988  and  Bq  >  5.5  kG.  The 

numerical  gridding  and  number  of  simulation  particles  is  the  sane  for 

the  two  runs.  The  result  indicates  a  much  better  choice  for 

initialization.  For  this  run  the  cutoff  radius  is  R  >  0.70  cm.  Taking 

the  scale  length  as  b  «  0.70  cm,  the  dimensionless  »  >  -0.11, 

2  2 

dimensionless  -  >27  and  1  -  -  1/Yq  >■  >012  .  The  maximum 

values  assumed  by  the  potentials  were  4/Yq  >  .069,  a^/rQ  «  .068, 

Sq/yq  -  .0027  (self-potential),  b^BQ/2YQ  -  .088. 

12.  C0HCLUSI(»IS 

Ve  have  derived  an  approximate  equilibrium  for  an  intense  electron 
bean  propagating  in  a  pure  solenoidal  field.  This  equilibrium  is  a 
higher  order  approximation  to  a  true  Vlasov  equilibrium  than  the  K-V 
equilibrium.  It  includes  space  charge  depression  effects,  the  shear  in 
the  longitudinal  nonentun  and  the  self-induced  longitudinal  diamagnetic 
field.  Simulations  confirm  the  theoretical  results.  Contemporary  and 
proposed  experimental  beams  have  parameters  for  which  these  effects  may 
be  measurable,  even  though  experimental  beans  tend  to  have  rounded 
profiles. 
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Figure  Ic.  Normelized  RMS  emittance  c  /c 

-3  ** 

(ca)  as  a  function  of  Oq  (cm  )  for  Tq  ■  9. 
The  other  paraseters  are  the  external^ 
solenoidal  field  Bq  -  4.5  kG  (dashed),  5.5 
kG  (solid),  and  6.5  kG  (dotted)  and  ■ 
.97,  .988  and  .993.  ^ 
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